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ABSTRACT: We present results from an inelastic neutron scattering experiment 
on the candidate Heisenberg spin ladder vanadyl pyrophosphate, (VO)2P207. We find 
evidence for a spin- wave excitation gap of Egap = 3.7 ± 0.2 meV, at a band minimum near 
Q = 0.8 A~^. This is consistent with expectations for triplet spin waves in (VO)2P2C)7 
in the spin-ladder model, and is to our knowledge the first confirmation in nature of a 
Heisenberg antiferromagnetic spin ladder. 

PACS Indices: 75.10.Jm, 75.25. +z, 75.30.Ds, 75.50.Ee. 



I.Introduction 



Many neutron scattering studies of the magnetic excitation spectra of Heisenberg spin 
chains have been reported since Haldane first suggested that integer- and half-integer spins 
might lead to qualitatively different physics. ^'^ This conjecture, that only half- integer-spin 
chains are gapless, is now generally regarded as experimentally confirmed. 

To gain insight into the physics underlying the surprising one-dimensional result, one 
may study systems intermediate between spin chains and two-dimensional systems, which 
are gapless for all spins. The simplest of these is the "Heisenberg spin ladder", which 
has isotropic couplings Jy along the chains and Jj^ between them. Although this system 
with both antiferromagnetic and ferromagnetic^^"^*^ interchain couplings has been 
the subject of considerable recent theoretical interest, to our knowledge no examples have 
previously been confirmed in nature. 

Copper nitrate, Cu(N03)2 • IH2O, was an early ladder candidate, since a ladder su- 
perexchange pathway was evident, and the specific heat and susceptibility could be de- 
scribed by a ladder model. However, theoretical studies showed that ladders and alter- 
nating chains lead to very similar thermodynamics, ^■'^"^"^ and the issue was finally de- 
cided in favor of an alternating chain by proton magnetic resonance^^ and elastic neutron 
scattering^^ experiments. Although copper hydrazinium chloride^^''^'' is not a complete 
spin ladder, it apparently does form a closely related antiferromagnetic "sawtooth", in 
which Cu"'"^ and Cl~ ions occupy alternate ladder sites. 

Recently two other candidate ladder systems have been discussed in the literature, 
a Sr„_iCu^+i02n series^^'^^ (expected to consist of coupled ladders with a frustrated 
"trellis" coupling-'^-'^), a La4+4„Cu8+2nOl4-l-8n series"^^ (another coupled-ladder system), 
and the isolated-ladder antiferromagnet vanadyl pyrophosphate, "^^""^^ (VO)2P207, which 
is the subject of this study. 
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A ladder configuration of S = 1/2 ions is clearly evident in the crystal structure 
of (VO)2P207 35,37,39 (^ggg especially Fig.2 of Ref.[39]), with spacings of 3.19(1) A between 
rung ions and 3.864(2)A between chain ions. The susceptibility of this material is fitted 
surprisingly well by an alternating chain model"^^ with Ji = 11.3 meV and J2 = 7.9 meV, 
but this need not eliminate the ladder model if the thermodynamic similarities previously 
noted for Cu(N03)2 • IH2O apply in the (VO)2P207 parameter regime as well. A recent 
theoretical study^ confirmed this similarity and found that the (VO)2P207 susceptibility 
was also accurately described by a ladder model with nearly equal interactions, J|| = 
7.76 meV and = 7.82 meV. The singlet-triplet gap at J±/J\\ = 1.0077 is given by^^ 
Egap/J\\ = 0.50759(1), corresponding to Egap = 3.94 meV with these parameters. 

One may test the ladder model of (VO)2P207 directly in a neutron scattering ex- 
periment, through a determination of characteristic propctics of the spin-wave excitation 
spectrum. In the ladder model with the parameters quoted above the singlet-triplet gap 
was predicted^ to be ~ 3.9 meV at a momentum transfer of Qa = tt (equivalent to 
Q = 0.813 in (VO)2P207). This one-magnon band reaches a maximum energy of 

16 meV near 0.3 and 1.3 A""^, and a crossing spin-triplet two-magnon band gives a 
secondary gap of 7.9 meV at Qa = and 27r and a broad plateau at 17-18 meV centered 
on Qa — TT. (See Fig. 6 of Ref.[8].) Structure factor calculations on finite clusters^^ indicate 
that the strongest neutron scattering should occur near Qa = tt, and higher levels above 
these two triplet bands should not be strongly excited. In contrast, a fit of the alternating 
chain model to the susceptibility predicts a somewhat higher gap of 4.9 meV, and the asso- 
ciated momentum transfer is problematical (it depends on an unspecified alternating-chain 
pathway) and need not have any simple relation to the Q — 0.813 A^^ expected in the 
ladder model. 
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II. Experimental Procedure and Results 



At present, single crystals of vanadyl pyrophosphate of a sufficient volume for neutron 
inelastic scattering experiments are not available, however several experiments have demon- 
strated that neutron scattering experiments on polycrystalline samples of one-dimensional 
magnetic systems can yield considerable information.^^ This is particularly true when us- 
ing time-of-ffight spectrometers, which can access a wide range of {Q, u>) simultaneously. 
The powder average of the correlation function observed experimentally is given by^^ 



5(Q,a;)=T(a;)|F(Q2)| 1 / ^ S{qu,q^,u;) dQ (1) 

where T{(jj) is the temperature factor 

T{uj) = (l - eM-r^/kBT)) , (2) 

F{Q) is the ionic form factor, and ^| and are the parallel and perpendicular projections 
of the total momentum transfer Q relative to the intrinsic ID magnetic axis. Clearly for 
any given Q all values of < Q will contribute to the scattered spectrum. Thus a mode 
which requires |gj|| = qi along the ID axis to be excited will first appear experimentally at 
Q = qii but will persist at values Q > qi, because these can excite systems with oblique 
orientations such that |^|| = qi. 

Measurements of the excitation spectrum of a polycrystalline sample of (VO)2P207 
were performed on HET, a direct geometry chopper spectrometer on the ISIS pulsed neu- 
tron source at the Rutherford Appleton Laboratory. A rotating Fermi chopper phased to 
the neutron pulse monochromates the incoming neutron beam. Banks of detectors at 4 
m and 2.5 m cover the angular ranges of (f) = 3° to 7" and ^ = 9" to 29" respectively. A 
further two banks at these lengths are situated at average angles of 118° and 130°, and 
are used for estimating non-magnetic background scattering. 
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The sample was prepared by Brody, in a manner similar to that described by Johnston, 
Johnson, Goshorn and Jacobson."^^ The final annealing was performed under flowing He at 
973 K for 7 days. For this experiment it was mounted on the cold finger of a closed cycle 
refrigerator in a sachet of aluminium foil. All data described hereafter were collected while 
the sample was at a temperature of 13 K, using incident energies of 100 meV, 35 meV, 25 
meV and 15 meV. 

Data collected simultaneously across the full angular range of the forward scattering 
banks using an incident energy of 25 meV were combined to produce the contour plot of 
the scattering intensity shown in Fig.l. The dispersion relation predicted by Barnes and 
Riera^ has been superimposed on the plot to compare the ladder model to the experimental 
data, and to demonstrate the effects of powder averaging, which spreads the scattering 
over larger Q values than the minimum required for excitation. The predicted cut-off 
of the scattering at approximately 18 meV energy transfer is clearly seen, and at lower 
energies there is evidence of scattering consistent with the theoretical dispersion relation, 
particularly in the region Q ^ 0.8 A~^; the one-magnon band can be seen as it passes 
through a minimum energy near Q = 0.8 A"-*^, where there is an increase in intensity as 
would be expected from both a peak in the density of states and in S{Q,lij). 

Using a lower incident energy of 15 meV the resolution is sufficient to resolve a gap 
between the elastic peak and the inelastic magnetic scattering. Fig. 2 shows data summed 
over an angular range oi (p = 14. 3*^ to 18.7'', which at an energy transfer of 3.7 meV 
represents a Q range of 0.72 A"-*^ to 0.88 A"-*^ (see Fig. 2 insert). The data have been 
corrected for kf/k^ and detector efficiency. The figure clearly shows the onset of inelastic 
magnetic scattering above an energy transfer of approximately 3 meV. The low-energy 
onset of the inelastic magnetic scattering has been fitted with a Lorentzian convoluted 
with the resolution function, yielding a value of 

Egap = 3.7 ± 0.2 meV (3) 
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for the gap. This is illustrated in Fig. 2 by a dashed line in the upper graph and by a solid 
line in the lower graph, which shows the data after subtraction of the elastic and quasi- 
elastic components. The smaller dashed line on the lower graph represents the resolution 
function, and suggests that the spectrum has an intrinsic width, as expected if we are 
exciting a band rather than a discrete level. 

The open symbols on the upper graph represent the data which was collected in the 
2.5 m detectors at (p — 130''. At these high angles, and thus high Q values, there is no 
magnetic scattering due to the fall-off of the ionic form factor, and the inelastic scattering 
arises from single and multiple phonon scattering processes. This data has been multiplied 
by an empirical scale factor of 0.2 (determined in previous experiments) to provide an 
indication of the amount of non-magnetic background scattering underlying the magnetic 
scattering data collected at low Q values. The non-magnetic inelastic scattering contribu- 
tion estimated in this way is clearly very small, and there are no features coincident with 
the observed magnetic scattering. The non-magnetic background scattering estimated in 
this manner has not been subtracted from the data presented in the lower plot. This 
data shows a feature at 1.8 meV energy transfer which is also visible in the scaled high 
angle data, confirming that the non-magnetic background scattering is indeed as weak as 
estimated by the scaled high angle data. 

The resolution function of the spectrometer is made up of a contribution from the mod- 
erator and a contribution from the transmission of the chopper. The former is described 
by a Gaussian convoluted with an exponential, the later simply by a Gaussian. The data 
below an energy transfer of 4 meV has been fitted using the FRILLS least squares fitting 
routine. The elastic line is well described by the resolution function, but there appears 
to be a second broad component under the elastic line, which has been fitted with a 
Lorentzian convoluted with the resolution function. The origin of this scattering is not yet 
understood, and may involve scattering from magnetic defects in the sample. 
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No attempt was made to fit the scattering intensity versus energy above the gap using 
the ladder model, because there is an obvious discrepancy; the observed scattering is 
considerably stronger than expected from numerical studies of S{Q, u) on small lattices, ^-"^ 
which anticipate a rapid fall in intensity away from the gap. This is the single feature of 
our data which does not appear to be in accord with current expectations for the ladder 
model. 

The overlap of the one- and two-magnon bands would give rise to considerable structure 
in the scattering above the gap energy, notably contributions from two-magnon production 
at a; 8 and 17 meV, for which we do see some evidence. This includes scattering above 
10 meV near Q = 0.8 which is evident in our 35 and 100 meV data. The 100 meV 
data indicates that this higher energy, Q ^ 0.8 A"-*^ scattering is strongest around 15 
meV, and may therefore be due to the two-magnon band expected at 17 meV. Due to 
the gaps in experimental coverage we cannot be certain about the interpretation of these 
higher-energy excitations. 

If data were available across a broader Q range, and without the null region evident in 
Fig.l, which is a consequence of the space between the 2.5 m and 4 m banks of detectors, 
it should be possible to extract the shape of the two branches and to study the indications 
of anomalous behavior in S{Q,U!). Such an experiment is planned for the near future. 

III. Summary and Conclusions 

In this note we present results from an inelastic neutron scattering experiment on the 
antiferromagnet vanadyl pyrophosphate. We find evidence for a band of triplet spin waves, 
with a minimum excitation energy of Egap = 3.7 ± 0.2 meV at a momentum transfer near 
Q = 0.8 A"-*^. These values are consistent with the expectations of the spin ladder model of 
(VO)2P207, which predicted Egap — 3.9 meV at Q = 0.813 A~^. The dispersion relation 
u;{Q) of spin excitations from approximately 0.6 A~^ to 1.5 A~^ is evident in our data, 
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and is also consistent with expectations for the fowest one-magnon band in the ladder 
model. These results confirm that (VO)2P2C)7 to a good approximation is a realization of 
the antiferromagnetic Heisenberg spin ladder. 

There remain two outstanding problems with the interpretation of neutron scattering 
from (VO)2P207. First, the intensity of scattering above the gap appears to be much 
stronger than expected from structure factor calculations on small lattices. Second, al- 
though the broad features of the excitation spectrum were observed in this experiment, 
the details of the expected bands above 10 meV are rather complicated, and an experimen- 
tal investigation will require more complete coverage in {Q,u). Ideally this would involve 
a scattering experiment using a single crystal. 
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Figure Captions 



Figure 1 

Contour plot of the scattering intensity from (VO)2P207 at 13 K using an incident energy 
of 25 meV. Contours of intensity run from 0.5 to 1.3 in steps of 0.1. (Arbitrary units.) 
The one-magnon (sohd hne) and two-magnon (dashed hne) bands predicted by Barnes and 
Riera^ are superimposed on the data. 

Figure 2 

A constant (j) scan across the predicted band minimum. Data from eight detectors which 
trace out the trajectories shown in the inset, have been summed. The upper graph shows 
the data (filled circles) and the fit as described in the text, together with data from the high 
angle bank of detectors (open circles), which have been multiplied by an empirical scale 
factor of 0.2 to provide an estimate of the non-magnetic background scattering. The lower 
graph represents the magnetic scattering after subtraction of the elastic and quasi-elastic 
components. 
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